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ABSTRACT 



Aims. H.E.S.S. observes an increasing number of large extended sources. A new technique based on the structure of the 
sky map is developed to account for these additional structures compared to the common point source analysis based 
on Li and Ma. 

Methods. Minkowski functionals are powerful measures from integral geometry. They can be used to quantify the 
structure of the counts map, which is then compared to the expected structure of a pure Poisson background. Gamma- 
ray sources lead to significant deviations from the expected background structure. The likelihood ratio method of Li and 
Ma is exclusively based on the number of excess counts and discards all further structure information of large extended 
sources. The morphometric data analysis incorporates this additional geometric information in an unbiased analysis, 
i.e. without the need of any prior knowledge about the source. 

Results. We have successfully applied our method to data of the H.E.S.S. experiment. The morphometric analysis 
presented here is especially designed for the detection of faint extended sources. 
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1. Introduction 



The aim of this study is to introduce a novel approach to 
data analysis in very high energy (VHE) gamma-ray as- 
tronomy, where extended sources are detected via morpho- 
metric valuations. 

Early studies in VHE gamma-ray astronomy focused on 
the study of point sources. For this purpose highly efficient 
analysis techniques were established, like the quantification 
of the significance of a photon count excess using a likeli- 
hood ratio method ( Li & Ma 1983 ) . But the increasing num- 
ber of large extended s ources (e.g 



Aharonian et al. 



2006a 



Aharonian et al.| 2007 ) and first d etection of diffuse VHE 
emissions ([Aharonian et al.||2006b) emphasize the need for 



new approaches, which might be more suitable for extended 
structures. 

An analysis based exclusively on the number of excess 
counts above the expected background level discards all 
information about the shape of the region where the ex- 
cess is observed. Furthermore, it does not use the informa- 
tion of possible correlations of nearby excess regions. While 
this additional information is negligible for point sources, it 
might provide a means to detect and study faint extended 
sources, which are too weak to be seen when looking at the 
amount of excess photons only. 

* e-mail: daiiiel.goering@physik.uni-erlangen.de 
** e-mail: michael.klatt@physik.uiii-erlaiigen.de 



A well known approach to include all available infor- 
mation in an analysis is the full likelihood fit of a model 
to the measured data, as used by high ene rgy gamma-ray 
i ^elescopes like EGRE T (I Mattox et al.|1996[ ) or Fermi/LAT 
(Atwood et al. 2009 1 . While likelihood analyses are very 
powerful, they require lots of a priori knowledge to build 
a proper model for the background and potential sources. 
The quality of a likelihood analysis is strongly influenced 
by the quality of the chosen models. 



This paper shows how to incorporate the additional 
structure information of extended sources into an analy- 
sis without the need for prior knowledge about the source. 
To achieve this, a reliable and powerful technique to quan- 
tify structures of gamma-ray counts maps is needed. The 
Minkowski functionals provide this techniq ue. They are well 
studied tools from integral geometry ( [Schneider fc Weil 
2008 Santal6|1976| ) and they are vastly used for both study 



ing structures in statistical physic s ( Mcckc 1998 Mecke & 
ISt oyan|[2000} !Schr6der-Tu rk et al.[[2010M20lTraiid for pat- 



tern analysis ( Mcckc 1996, [Becker et al |2003 Mantz et al 
200^. T hey were already successfully applied in astron 



omy (e.g. Mecke et al.|1994 ). They were used to investigate 
point processes in cosmology and the large scale structure 
of the univer se (Kcrscher et al.[[2001 Kerscher et al.[[2001 



Colombi et al._ 2000J and as probe s of non-Gaussiamty m 



the cosmic microwave background (Schmalzing et aL][1999 



Gay et al.[[2012{[Ducout et al.[[2013 ) 
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Fig. 1. Structure quantification via Minkowski function- 
als. (a) Counts map, simulated Poisson distributed random 
number of counts k. In order to characterize the morphol- 
ogy, the image is turned into a black and white image via 
thresholding — see (b); the three Minkowski functionals are 
then evaluated for the b/w image, (b) Area A (c) Perimeter 
P (d) Euler characteristic x 



This work uses for the first time Minkowski function- 
als to detect sources in gamma-ray astronomy by quanti- 
fying structures in gamma-ray counts maps. The paper is 
organized as follows: section [2] introduces the Minkowski 
functionals in detail and describes the quantification of the 
structure of the counts map. The morphometric analysis, 
i.e. the actual source detection via structure quantification, 
is presented in section [31 How to express the background 
morphology and how to define a compatibility of the struc- 
ture of a measured counts map with the background is ex- 
plained in section |3.1| using a global null hypothesis. The 
technique is extended in section [3. 2| to local structure devi- 
ations found with Minkowski sky maps allowing to resolve 
and localize the gamma-ray sources. Section [3 .3| applies the 
analysis to simulated data. Finally, the results for counts 
maps observed with H.E.S.S. are given in section [i] 

2. Structure characterization 

This section shows how to characterize the structure of a 
gamma-ray counts map. While similar methods may be 
used to quantify the morphology of extended gamma-ray 
sources, this is not the subject of this paper. Although 
the following structure analysis has not yet been applied 
in gamma-ray astronomy, it is well studied in integral ge- 
ometry and in statistical physics dSchneider fc Weil 2008 



Mecke|p98{ [Mecke fc Stoya^|2000 

A gray-scale image, here the counts map, is turned into 
a black and white (b/w) image (Mecke 1996) . For each 
threshold value p, all pixels with counts k > p are set to 
black, the others remain white — see Fig. [l] The structure 



Table 1. Look-up table for Minkowski functionals: the 
functional values of area A, perimeter P and Euler charac- 
teristic X are assigned to each 2x2 neighborhood of the im- 
age. The unit of length is t he edge-length of a pixel. Simila r 
data can be found in Refs. ( |Mecke|1996{|Mantz et al.|2008[ ). 
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of the image is then analyzed as a function of the threshold 

P- 

The structure of each b/w image is quantified by the 
Minkowski functionalqil In two dimensions there are three 
of them. They are proportional to well known geometric 
quantities: the area A of the black pixels, their perimeter P 
and the Euler characteristic x> which is the integral of the 
Gaussian curvature. It is a topological constant; for closed 
domains it is given by the number of components minus the 
number of holes. Figure [l] vi suali zes how a counts map (a) 
is turned into a b/w image \{h]' " 



which is then quantified 
The area as a function 



IS turned into a b/w ima ge \(b, 
by Minkowski functionals |(b)|(( 



of the threshold contains the knowledge about the num- 
ber of counts. However, it does not supply any information 
about their arrangement, for which additional information 
is provided by the perimeter and the Euler characteristic. 
The Minkowski functionals are powerful shape mea- 
sures. Due to their additivity and continuity, they are ro- 
bust against noise and have short computation times. There 
are several linear time algorithms for the calculation of area, 



perimete r and Euler chara cteristic (e.g. Mantz et al.||2008 
Schrode r- Turk et al.|2010 ) and for 3d data (e.g. Arns et aL 
|2010[|ScEroder-Turk etal.|2011| [Schroder-Turk etal.|2012p 



The straightforward algorithm used here is based on the 
look-up table [l] The image is decomposed into 2x2 neigh- 
borhoods. To each of the 16 possible configurations the 
value of the Minkowski functionals is assigned. Due to their 
additivity the sum of the local contributions yields their 
global value. The unit of length is defined as the edge-length 
of a single pixel, thus the area of a pixel is one. In order 
to avoid multiple countings when iterating over the whole 
image, only the part may contribute which is unique to a 
2x2 neighborhood, i.e. each quarter of the four pixels next 
to the center. For example, a single black pixel has area 
and Euler characteristic one and perimeter four. However, 
when iterating over the image, it will appear in four differ- 
ent 2x2 neighborhoods, namely configurations two, three, 
five, and nine. Thus, Table [T] assigns to each of them area 
and Euler characteristic one fourth and perimeter one. The 
white pixel in configuration 15 can be interpreted as part 
of a hole; it contributes negatively to the Euler character- 
istic. In configurations seven and ten in Table [T] the black 
pixels sharing only a vertex are chosen to be connected. If 



^ Other names are valuations, Quermafiintegrals, intrinsic vol- 
umes or Hadwiger measures. 
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they were disconnected, the weights for the Euler charac- 
teristic would be positive. The choice is arbitrary, as long as 
the probability distribution for the Euler characteristic is 
calculated consistently. However, connecting them helps to 
distinguish a single cluster of black pixels from two domains 
distant to each otherQ 

The choice of boundary conditions has a strong influence 
on the stru cture quantification and its efficiency [Stoyanj 
et al. 1987). Throughout this work, closed boundary con- 
ditions are applied. This means all pixels outside the win- 
dow of observation are set to white, thus all domains are 
closed. A discussion of the different impacts and drawbacks 
of various boundary conditions goes beyond the scope of 
this paper and will follow in future publications. 

Hadwiger's completeness theorem ensures that the 
Minkowski functionals provide a robust and comprehensive 
morphology analysis, i.e. the Minkowski functionals form a 
complete basis of all valuations which are defined on unions 
of convex sets and which are motio n invariant, addit ive and 
at least continuous on convex sets (IHadwiger"1957'). In or- 
der to quantify anisotropy, they can be generalized to tensor 



2010 2011). 



valuations (e.g. Schneider & Weil 2008 Schroder- Turk et al. 



3. Source detection 

3.1. Global null hypothesis test 

Structural deviations from the background morphology 
are to be detected. Therefore, the characteristic structure 
of a background measurement has to be known. As the 
Minkowski functionals quantify the morphology, the prob- 
ability distribution P of their functional value X e {A, P,x} 
for a counts map of a background measurement is needed. 
The first step is to choose a suitable background model. 

In general, determining fiX) is the most complex task in 
the morphometric analysis. A reasonable model for back- 
ground counts in VHE gamma-ray astronomy is the as- 
sumption of homogeneously and isotropically Poisson dis- 
tributed counts in each bin of a sky map with equal area 
bins. This is due to the fact that the largest part of back- 
ground events in ground based VHE gamma-ray astronomy 
is caused by VHE hadrons. These hadrons loose their di- 
rection correlations in interstellar magnetic fields and ar- 
rive at the earth as a uniform flux of VHE particles from 
every direction. For a real measurement the homogeneous 
and isotropic background will of course be distorted by de- 
tector effects and non-uniform exposure of the sky. As will 
be shown later in section |4.1[ the data can be corrected 
for these effects, i.e. the typical structure of a background 
measurement can be derived from the structure of a pure 
homogeneous and isotropic Poisson backgroundr] 



^ If a marching square algorithm is used to find a more com- 
plex triangulation of the domain of black pixels, the weights 



for ar e a and perimeter have to be adjusted - see Ref. |Mantz| 
(|2008L The probability distributions for the Minkowski 



et al. 



functionals have to be calculated consistently. However, no sig- 
nificant effect on the final results has yet been observed. 

'^ The choice of bin size is free, because a uniform Poisson field 
remains a homogeneous Poisson field for any chosen bin size and 
an arbitrary point spread function; the null hypothesis will be 
unchanged. However, to quantify the actual source morphology 
and not random noise, the bin size should be adjusted to the 
point spread function. 
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Fig. 2. Structure of background noise for a 5 x 5 bin field 
at different thresholds p resulting in different probabilities 
Pp - YjT^p ^' ^ '*/'! for a black pixel. Probability distributions 
V of (a) area A; (b) perimeter P] (c) Euler characteristic ;if. 
(d) Compatibility C of Euler characteristic ;^f, with Cix - 
-1) the sum of all probabilities P < P(x = -1), which is 
visualized by the equally large colored regions beneath the 
graphs in (c) and (d) - mind the different scales. 



The likelihood for the area A of the black pixels is a 
suitable introductory example; the probability distribution 
for each threshold p is given by the Binomial distribution 

PiA) = Bmom(A,pp,N^), 

where k is the number of black pixels, A^^ the number of pix- 
els in a A^ x A^ sky map and pp the probability that a pixel 
is black, i.e. the counts are larger or equal than the thresh- 
old p. Assuming a Poisson background with an expected 
number of counts per bin of A, pp is given by 



Pp 



z 



(1) 



The probability distributions 'P(A), P{P), and P{x) of the 
area, perimeter, and Euler characteristic, respectively, are 
plotted in Fig. [2] (a-c) as a function of the probability pp 
that a pixel is above the threshold p. For small A^ these 
distributions can be found by evaluating all possible b/w 
images and inferring the distribution by counting equivalent 
pixel configurations. 

Gamma-ray sources can be detected by looking for 
structures, which are very unlikely, given the hypothesis 
of a pure background measurement was true, i.e. assuming 
that there are only background events within the observa- 
tion window. A probability measure that a given counts 
map with a structural value X is compatible with this hy- 
pothesis may be defined as follows: The compatibility 



CiX) = Yj ^(^') 

'P(Xi)<'P(X) 



(2) 
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is the probability that a less likely structure appears. 
Figure [2] (d) shows the compatibility C(x) for the Euler 
characteristic from Fig. [2] (c) . 

The definition of the compatibility follows the scheme 



10 



given m 
efficient 



Neyman & Pearson (1933) to construct a most 
lypothesis test. The form used here may be de- 
rived from the general scheme given in the paper, by setting 
the supremum of alternative hypotheses to 1, i.e. impos- 
ing no constraints on alternative hypotheses. The hypoth- 
esis of a pure background measurement shall be rejected, 
if the compatibility is less than 0.6 ■ 10"*'. This hypothesis 
criterion is adjusted to the commonly used 5cr deviation; 
a normally distributed random variable deviates from the 
expected value at least 5 cr with a probability of approxi- 
mately 0.6 -10"^ 

Instead of dealing with tiny compatibilities it is often 
more convenient to use the logarithm of this likelihood value 
or to define the deviation strength 



D:^- logio(C) 



(3) 



The conversion between compatibility C and standard de- 
viation cr is given by Cicr) - erfc (cr/ V2], where erfc(.j;) = 
^ r exp(-f^) df is the error function. 

Figure [3] depicts a structure analysis both of a rather ex- 
tended and of a more point like source, investigating simu- 
lated data of Poisson distributed random number of counts. 
The deviation strength D is plotted over the threshold p. 
As for both sources there are thresholds for which the de- 
viation strength is greater than 6.2, the null hypothesis of 
a pure background measurement can be rejected in either 
case. The results differ for the two sources, while their total 
flux is equal, and reveals basic information about the ex- 
tension of the sources. The point source cannot be detected 
at thresholds near the background level A and is only ap- 
parent for high thresholds, i.e. only bins with high counts 
contribute to the detection. For the extended source only 
counts in the order of the magnitude of the background 
fluctuations are present in the sky map. However, as the 
number of bins with slightly increased counts exceeds the 
typical background predictions, the hypothesis is rejected 
nevertheless. 

3.2. Local Minkowski sky maps 

So far, it is possible to scan the entire field of view (FoV) 
of an observation for additional structures on top of the ex- 
pected background. In order to localize the deviations and 
thus, locate gamma-ray sources and gain insight in their 
extension and morphology, further information is needed. 
Due to their motion invariance, the scalar Minkowski func- 
tionals can not be used directly to localize structuresjj 

Instead of analyzing the entire FoV with the methods 
introduced so far, a small sliding window may be used. 
This window can be moved across the FoV to study the 
local structure of the sky map in the sliding window. Using 
this approach it is possible to construct Minkowski sky 
maps from the counts maps, which are just as useful as 
significance maps constructed using the approach from[Li| 



■* While the generalization to motion covariant tensorial val- 
uations might provide this insight, these generalization will not 



be discussed in this paper. Instead a straightforward and easy 
approach to the problem is discussed here 
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Fig. 3. Deviation strength D{A) for structure quantification 
via area A as a function of the threshold p for an extended 
and for a point source of equal flux within an observation 
window of 100 x 100 bins. 



& Ma (19831, but provide the additional sensitivity from 
the structure information used to determine the deviation 
strength. 

Since the maximum deviation strength is assigned to a 
pixel, a trial factor must be added to the result in order to 
avoid an overestima tion of the significance of found devia- 
tions (Goring 2012). Each of the different b/w images for 
different thresholds p contributes a trial to the search for 
structure deviations and as the number of trials increases, 
the probability to find a significant random fiuctuation in- 
creases. Suppose you look for measurements with a com- 
patibility of less than a, then the probability to find such 
a deviation in one image is a. The probability to find no 
such deviation in t independent images is (1 — a)' and the 
probability to find at least one such deviation in the image 
set is 1 - (1 - a)'. Thus, the significance after t independent 
trials is linked to the pre-trial significance a via 



l-(l-ay ^ta + 0(a^). 



(4) 



So, for a <K 1 the influence of t independent trials may be 
approximated by the so called trial factor, i.e. by multiply- 
ing a with t. For the corresponding deviation strength D, 
this results in an offset of logio(f) compared to the pre-trial 
deviation strength D, i.e. WD » : £), ~ £) - logjoCf). 

Obviously the different b/w images resulting from 
thresholding are not statistically independent, as they orig- 
inate from the same intensity profile. If the structure of 
a b/w image at a certain threshold can be interpolated 
from the structures of the b/w images of the next higher 
and lower thresholds, it does not contribute a separate 
trial to the search for deviations from the null hypothe- 
sis. Therefore, if t is set to the number of threshold steps, 
Eq. ^ will yield a conservative estimate of the post-trial 
significance. 

A rough estimate of the systematic error on the devi- 
ation strength D introduced by ignoring the infiuence of 
trial factors is given by 21ogi0(A^). This estimate is based 
on the fact that the number of different b/w images after 
thresholding is equal to the number of different gray levels 
in the original gray scale image representing the gamma- 
ray counts map. There are at most A^^ different gray levels 
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Fig. 4. Computation of a Minkowski sky map. LEFT: a 
given counts map. MIDDLE: for each pixel the local struc- 
ture within the sliding window is characterized (top) and 
the maximum deviation strength for all thresholds is as- 
signed to it in the Minkowski sky map (bottom) . RIGHT: it- 
erating through all pixels provides the Minkowski sky map. 



in an image oi N x N pixels and thus at most A^^ different 
trials may contribute to the search for structure deviations; 
from t < N^ follows 



D, *£)-logi„(r)>£)-21ogio(A^). 



(5) 



Therefore, £)-21ogig(A^) may be used to compute a conser- 
vative estimate of the post-trial deviation strength; for the 
preceding discussion see Ref. ( G6ring|["2012 ). 

Figure |4] depicts the concept of a Minkowski sky map. 
For each of the inner pixels of the counts map the lo- 
cal structure is to be characterized (depicted on the left). 
Therefore, the Minkowski functionals of a sliding window 
with A^xA^ pixels are evaluated (illustrated in the top picture 
in the middle of Fig. El) . The maximum deviation strength 
for all thresholds p is assigned to the pixel at the center of 
the sliding window (plotted at the bottom picture in the 
middle). Iterating over all inner pixels for which the sliding 
window is completely within the counts map provides the 
Minkowski sky map. 

The information content and intuitive interpretation of 
such Minkowski sky maps can be enhanced, by adding a 
sign to the deviation strength of the different pixels. By 
choosing the sign to be negative if A < N^p, i.e. if the num- 
ber of black pixels is smaller than expected, and positive 
otherwise, the sign of a pixel shows if the local deviation is 
due to an overestimation of the background or due to ad- 
ditional flux from potential gamma-ray sources. Minkowski 
sky maps detect local structural deviations and depict them 
in a both illustrative and quantitative image, depicting the 
lack of trust in the hypothesis that there are only back- 
ground fluctuations. 

Although the null hypothesis is tested locally, the back- 
ground assumption is a global null hypothesis, i.e. the ex- 
pected number of counts per bin A is chosen globally. If A 
is chosen locally for every sliding window, structures larger 
than the sliding window may result in an increased A and 
no deviation at all if the local structure is homogeneous and 
isotropic within the sliding window. 




Fig. 5. Structure analys is o f a simulated sky maps 
Given intensity proflle. (b) Simulated counts map. 



^ 



Signiflcance map, evaluate d w ith standard techniques 



weighting the count excess, (d) Minkows ki sk y map, the 
structure is c haracterized with the area A, |(e)| the perime- 
ter P and 1(f) I the Euler characteristic ;^f. 



3.3. Simulated data 

It is now possible to detect, to localize, and to study the 
shape of gamma-ray sources using the morphometric anal- 
ysis with Minkowski sky maps, introduced in the previous 
section. Figure [S] depicts an analysis of exemplary simu- 
lated data. Figure (a) shows the given intensity profile of the 
test pattern sources, the expected number of counts per bin 
caused by the chosen sources. The test pattern consists of 
Gaussian shaped extended sources with a ratio of the semi 
major axes of 2. The peak intensities of the sources along 
a central ray from the image center are equal. The sizes of 
sources with the same distance from the center are equal 
and the peak intensities are increasing in counter clockwise 
direction. Plot (b) is the simulated counts map analyzed 



in (c -f). In figure (c) the standard analysis technique oflLi 
fc M a (1983 J as used by the H.E.S.S. experiment (cf. F 
Aharonian et al.||2006 ) is applied to the counts map (b) for 
comparison, l^'igures (d-f ) are the Minkowski sky maps with 
a sliding window size 5x5; the structure is quantified by 
either the area A, the perimeter P or the Euler character- 
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istic X- The same 5x5 window i s used as ort-reg ion for the 
significance determination using Li & Ma ( 1983 1 . 

The first apparent resuh drawn from Fig. [5] is, that the 
single functionals and the standard analysis are similarly 
sensitive to sources; area, perimeter and Euler characteris- 
tic are comparably competitive in finding gamma-ray sig- 
nals. 
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4. H.E.S.S. data 

4.1. Detector acceptance correction 

An analysis of real data must always take into account 
the detector acceptance. After modeling the camera accep- 
tance, the effect must be corrected in order to regain an 
isotropic and homogeneous structure for background mea- 
surements. For each bin / only a fraction /■ of the signals 
are expected to be detected. If each bin is simply weighted 
with l/fi, fractional photon counts will occur destroying the 
Poisson structure of the sky map. 

The null hypothesis that for an ideal camera acceptance 
there is only background noise with intensity A, allows an 
acceptance correction, which preserves the Poisson struc- 
ture. As only the fraction /■ of the events in bin / are de- 
tected, the actual background intensity Aj - fi ■ A varies 
with each bin. Following the null hypothesis, that the num- 
ber of counts is a Poisson distributed random variable with 
mean /l,-, the original random process with mean A can be 
regained, if a new Poisson distributed random variable with 
mean A^ - {I - fi) ■ A \s added. 

In the center of the field of view, where /■ ~ 1, the num- 
ber of counts remains effectively unchanged. For /■ <K 1 the 
additionally created pseudo photon counts may cover fea- 
tures, but they never introduce additional structural devia- 
tions from the homogeneous isotropic Poisson field, as they 
fulfill the given null hypothesis by construction. Covering 
regions of low acceptance with a layer of pseudo events cor- 
responds to the fact that the instrument is less sensitive to 
signals in these regions compared to regions with high ac- 
ceptance - see also Refs. ( [Goring 2008 , Klatt 2010 Goring 
20121. 



4.2. H.E.S.S. source 

Using the methods discussed so far it is possible to an- 
alyze experimental data of ground based VHE gamma- 
ray telescopes like the H.E.S.S. experiment. As the mor- 
phometric analysis is targeted at extended structures, 
RX J1713.7-3946 is chosen for an exemplary analysis. This 
source has one of the largest angular diameters of the 
sources seen by H.E.S.S. and its morpho logy has already 
been studied in detail and is well known (Aharonian et al. 
2006a), which makes it a suitable benchmark for further 



analyses. 

The data set used for analys is corresponds to the data 
used in Aharonian et al. (2006a), but instead of the Hillas 
based event reconstruction discussed there the advanced 
event rec onstruction based on a likeli hood model fit pre- 
sented in de Naurois & RoUand ( 2009 ) is used to compute 
the list of reconstructed events from the recorded cam- 
era images. These reconstructed events are used to fill the 
binned counts map, which is the main input of the morpho- 
metric analysis. 



Fig. 6. Minkowski sky maps of RX J1713.7-3946. The 
morphometric analysis is based on (a) area A, (b) perimeter 
P and (c) Euler characteristic x- 



The second input needed is an acceptance map describ- 
ing the spatial sensitivity of the given observations, which is 
used to perform the previously discussed acceptance correc- 
tion. This map is created using standard H.E.S.S. analysis 
tools. In particular, the so called 2D acceptance model dis- 
cussed in|de Naurois|(|2012|) is used. The overall background 
level A is 



de Naurois 

determined from the counts map by excluding all 
regions containing known sources of VHE gamma-rays and 
computing the mean of the counts in the remaining bins 
normalized by the corresponding acceptance. 

Figure |6] shows the resulting Minkowski sky maps of a 
morphometric analysis based on the area A, the perime- 
ter P and the Euler characteristic x respectively. The un- 
derlying counts and acceptance maps used for this analysis 
were created using square bins of 0.02 ° width. The resulting 
sky maps clearly show RX J 17 13.7 -3946 and are in good 
agreeme nt with the results of the s tandard H.E.S.S. analysis 
given in Aharonian et al. (2006a). This demonstrates that 
an analysis based on the structure of the measured counts 
map is indeed possible and provides results comparable to 
those of well established tools. 



5. Conclusion and outlook 

The introduced morphometric analysis provides a novel ap- 
proach to data analysis in VHE gamma-ray astronomy. In 
contra s t to t he commonly applied hypothesis test by |Li] 
& Ma (1983) to detect a significant excess of gamma-ray 



counts on top of the expected background, it allows to 
incorporate additional morphometric information into the 
analysis of gamma-ray counts maps. Still, the underlying 
model depends only on A, the expected background level of 
a given measurement, and there is no need for a priori mod- 
eling of potential sources, as opposed to an analysis based 
on a likelihood fit of a comprehensive model to the data as 
commonly used in s atellite based gamm a-ray experiments 
like e .g. Fermi/LAT ( [Atwood et al.[2009| ) . Ref. ( |Klatt et al.' 
2012 | provides a short introduction to the shape analysis 
of counts maps. 

The given morphometric analysis is based on the char- 
acterization of the typical structure of a pure homoge- 
neously and isotropically Poisson-distributed background 
counts map. This typical structure is determined using the 
Minkowski functionals — morphometric valuations from 
integral geometry. Significant deviations from the typical 
background structure in measured gamma-ray counts maps 
can be used to detect gamma-ray sources in the same way 
as significant excess counts are used in analyses based on 
Li_&Ma (1983j). 

Using the presented basic ideas it is already possible to 
qualitatively reproduce analysis results based on well es- 
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tablished tool chains. In order to quantify the agreement 
of the different analysis techniques and potential sensitiv- 
ity gains of the new morphometric approach further in- 
depth studies are required. Still, there are strong indications 
that the presented methods will lead to a significant sen- 
sitivity gain in the foreseeable future. Already the current 
proof-of-concept studies show an impressive performance 
and a similar sensitivity of the three Minkowski functionals 
to deviations from the background structure. Combining 
the different functionals to refine the background charac- 
terization may well lead to a significant sensitivity boost. 
Furthermore, the Minkowski functionals can be generalized 
to tensor valued valuations. These Minkowski tensors quan- 
tify additional morphometric information like isotropy and 
homogeneity and thus provide a comprehensive view of the 
available morphological information. Incorporating this ad- 
ditional information may increase the sensitivity of the mor- 
phometric analysis even further. The inherent potential of 
the presented methods provides some exciting perspectives 
for data analysis in VHE gamma-ray astronomy. 
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